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Is the relatively isolated member of the neutrino mass spectrum heavier or lighter than the two 
closely-spaced members? This question — the character of the neutrino mass hierarchy — is of 
great theoretical interest. All previously identified experiments for addressing it via neutrino oscil¬ 
lations require that the currently unknown size of the U e 3 element of the leptonic mixing matrix 
(parameterized by the unknown #13 mixing angle) be sufficiently large, and will utterly fail in the 
limit #13 —> 0. For this reason, we explore alternative oscillation approaches that would still succeed 
even if #13 vanishes. We identify several alternatives that require neither a nonzero U e 3 1 nor the 
presence of significant matter effects. All include multiple percent-level neutrino oscillation mea¬ 
surements, usually involving muon-neutrino (or antineutrino) disappearance and very long baselines. 

We comment on the degree of promise that these alternative approaches show. 

I. INTRODUCTION 

Neutrino experiments have revealed that neutrinos have mass and that, similar to quarks, leptons mix. Neutrino 
masses were discovered with the observation that neutrinos change flavor as prescribed by the mechanism of mass- 
induced neutrino oscillations [jj. 

Three-flavor neutrino oscillations suffice to explain all neutrino data, with the exception of those from LSND . The 
LSND result, however, still awaits confirmation from another experiment and will not be considered here. Nonetheless, 
even under the conservative assumption that there are no other surprises in the neutrino sector, our understanding of 
leptons remains incomplete. While we currently know the absolute value of the two neutrino mass-squared differences 
and two out of four mixing parameters that govern neutrino<->neutrino and antineutrino<-s-antineutrino oscillations, 
we have little or no information regarding* 

1 . the magnitude of the C/ e 3 -element of the lepton mixing matrix, or the value of # 13 ; 

2. the value of the complex phase S which parameterizes CP- (and T-)invariance violation in neutrino<-s-neutrino 
and antineutrino<->antineutrino oscillations; 

3. the neutrino mass hierarchy, or the sign of Am 2 3 ; 

4. the magnitude of the neutrino masses (we only know mass-squared differences and various upper bounds on 
other combinations of the neutrino masses); 

5. whether neutrinos are Majorana or Dirac fermions. 

Items 1, 2, and 3 can be addressed by next-generation neutrino oscillation experiments. Here, we wish to concentrate 
on item number three and the outstanding question: how can one experimentally determine the neutrino mass 
hierarchy? 

The neutrino masses mi , i = 1,2,3, are defined in the following convenient way: m\ < m|, and 0 < Am 2 2 = 
777-2 — m i < |A7b 2 3 |, where Amf 3 = m§ —m\. A positive value of Am 2 3 implies mj > m 2 and a so-called normal mass 
hierarchy, while a negative value of A mf 3 implies to 2 < m\ and a so-called inverted mass hierarchy. Examples of an 
inverted and a normal neutrino mass hierarchy are depicted in Fig.^ 

The question of whether the neutrino mass hierarchy is normal or inverted is of considerable theoretical interest 4]. 
For example, one can argue on general grounds that grand unified theories favor a normal neutrino mass hierarchy |5). 
This may be understood qualitatively by noting that grand unified models relate leptons to quarks, and the quark 
mass hierarchies are normal. An inverted neutrino mass spectrum, on the other hand, would be quite un-quark like, 
and would probably reflect an underlying global leptonic symmetry (see i4j and references therein) that gives rise to 
the near degeneracy of mi and m 2 . 


We adopt the “standard” definition of neutrino mass eigenstates (discussed below) and the lepton mixing matrix U 



2 


normal hierarchy 


inverted hierarchy 



(m,) 


□ 


(m,r 


(m,r 



FIG. 1: Two distinct neutrino-mass hierarchies that fit all of the current neutrino data, from [lj. The color coding (shading) 
indicates the fraction |t/ a i | 2 of each distinct flavor v a , a = e, fj,,T contained in each mass eigenstate Vi. i = 1,2, 3. For example, 
|t/ e 2 1 2 is equal to the fraction of the {m 2) 2 “bar” that is painted red (shading labeled as V e ’). 


It is easy to understand, qualitatively, why the current data are unable to distinguish positive from negative Ato 2 3 . 
Almost all information we have on A m\ 3 comes from disappearance at atmospheric experiments and K2K. These 
provide information on the survival probability P ^ which is given (in vacuum) by 

P^~ 1 - sin 2 26 » 2 3 sin 2 > ( L1 ) 

modulo terms proportional to A m\ 2 or \U e3 \. These are suppressed because A m\ 2 ^ |A. 77 T.f 3 1 (at the three sigma 
level, Am 2 2 /| Ato 2 3 | < 0.06 Q), which we know from solar and KamLAND data, and |C/ e 3 1 2 <C 1 (at the three sigma 
level, \U e3 \ 2 < 0.051 @), which we know from reactor and atmospheric data. Current experiments are not sensitive 
to these sub-leading terms, and it is clear that Eq. CD is only sensitive to Aro'f 3 |. Hence both a normal and an 
inverted mass hierarchy fit the data equally well. 

In this paper, we explore several different oscillation probes of the neutrino mass hierarchy, and discuss under what 
conditions these can be experimentally accessible. More specifically, we concentrate our efforts on probes other than 
“matter-affected 13 oscillations,” which have been heavily explored in the literature and are defined and very briefly 
discussed in Sec. El To date, all phenomenological studies of how to determine the neutrino mass hierarchy via 
neutrino oscillations have, to the best of our knowledge, relied on a non-zero value for | C/ e 3 1 , and almost all appeal to 
matter effects in the 13 sector. One of the goals of this paper is to argue that neither matter effects nor a nonzero 
value of \U e3 \ are required in order to determine the neutrino mass hierarchy via neutrino oscillations — at least in 
principle. 

There are many reasons for pursuing novel oscillation probes of the neutrino mass hierarchy. One is that we would, 
ultimately, like to compare distinct determinations of the mass hierarchy. Such a comparison provides a test of 
whether there is new physics in the lepton sector besides three flavor oscillations. Another (and, perhaps, the most 
important) reason for pursuing different hierarchy-exploring channels is that, if 0\ 3 turns out to be too small, the 
traditional approach will not work, so new experimental probes are not only desirable but absolutely necessary! 

Our presentation is the following. In Sec. El we briefly review how the mass hierarchy can be revealed if \U e3 \ is 
large enough and there are nontrivial matter effects at work. In Sec. EH we remind readers that even in the case 
of pure vacuum oscillations the mass hierarchy is, in principle, observable. We then review one realistic method for 
determining the mass hierarchy using only vacuum disappearance channels. In Sec. m we study in detail the case 
013 = 0, and discuss various approaches one might consider in this case to determine the mass hierarchy. In Sec. Ivl 
we summarize our results and offer some concluding thoughts. A technical presentation of the oscillation probabilities 
in constant matter for vanishing $13 is included as an appendix. 

Our approach is “theoretical,” in the sense that we study the behavior of oscillation probabilities as a function of 
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the neutrino energy E and the baseline L instead of the expected number of observed events for a given neutrino 
beam, detector type, size, resolution, etc. Nonetheless, our approach allows one to estimate “how precise” certain 
measurements need to be so that the inverted and normal mass hierarchy lead to experimentally distinct results. 


II. UNVEILING THE MASS HIERARCHY THROUGH MATTER EFFECTS IN 13 OSCILLATIONS 


One very promising way of revealing the character of the neutrino mass hierarchy is to probe oscillations which 
are significantly affected by matter effects in the 13 sector nil In this section, we very briefly describe one 
manifestation of this well known phenomenon. Readers are referred to nil for more details. 

By “matter affected 13 oscillations” we mean neutrino oscillations that take place under the following two conditions: 
(i) | U e z 1 2 is large enough* and (ii) the propagating medium contains electron number densities N e such that the matter 
potential, \Z2GfN b , where Gf is the Fermi constant, is comparable to the 13 frequency, \/Am\ 3 \/2E. Note that, while 
the latter can be arranged (at least in principle) by choosing an appropriate neutrino beam and propagating medium, 
the former requires cooperation from Mother Nature. 

Matter affected 13 oscillations can be studied, for example, in accelerator-based, long-baseline searches for —> v e 

0. In constant matter densities, the —> v e oscillation probability P^ e is given by 


Pfie - PeF - sin 2 0 23 sin 2 20 ff sin 2 , 


e/i. 


'13 


Q 4 n 2 oneff _ A- 3 sin 2 d 13 
Sin ZU 13 — (A eff)2 ! 

A^f = yj (Ai 3 cos 26*13 - A ) 2 + Af 3 sin 2 20 i3 , 

A13 = Am * 3 
^ 13 2 E ' 


(II.1) 

(II-2) 
(II.3) 
(II-4) 


where A = \Z2GpN e is the matter potential. Eq. ITT. 1 1) is a good approximation as long as A 12 terms can be ignored. 
A®f is sensitive to the relative sign between Amf 3 cos 20 13 and A. Since A, cos 20 13 > 0, a normal neutrino mass 
hierarchy implies that, for a fixed value of $ 13 , P™ ax is larger than its vacuum counterpart, sin 2 O 23 sin 2 26 * 13 , while 
Pfie is smaller .' The opposite behavior is expected for an inverted hierarchy. 

While P Me is sensitive to the mass hierarchy, measuring it alone at approximately fixed values of L and E is often not 
enough to experimentally determine the mass hierarchy. This is due to a well known 0 i 3 -mass hierarchy degeneracy, 
which can be understood in the following way. If one only measures sin 2 20® 3 , for every normal hierarchy characterized 
by (Amf 3 , 6 *i 3 ) = (Ara^ , 0^), there will be an inverted one with (Amf 3 , 6 * 13 ) = (A mf 3 , 6 33 ) that yields exactly the 
same value of sin 2 20ff. Note that, as long as Ato 2 ^" ~ -Am^") the degeneracy is dominated by the relative size of 
0and 0 33 and the sign of Am 2 3 . Qualitatively, sin 2 20ff can be chosen the same for a normal and an inverted mass 
hierarchy if one chooses 9 + 3 to be less than 0 33 . 

The degeneracy described above can be broken by more experimental information. Simple ways of determining 
the mass hierarchy include comparing P^e and P^ g at A 33 P ~ 7r (which can be done with one detector and both 
muon-neutrino and muon-antineutrino beams), or comparing measurements of P Me at A^L ~ 7 r for different values of 
L (which requires two different detectors, but only muon-neutrino beams). Another way is to independently measure 
sin 2 20i 3 precisely enough such that one can statistically distinguish between 013 - 013 and 0 i 3 = 0 13 . 

Other probes of the neutrino mass hierarchy that depend on matter effects and a nonzero \U e z\ include precise 
atmospheric neutrino studies Q and measurements of the neutrino flux emitted during a nearby supernova explosion 

a 


The procedures referred to above do not take advantage of the fact that, in the case of three neutrino flavors, neutrino 
oscillations are governed by two distinct oscillation frequencies. Indeed, extracting the neutrino mass hierarchy via 
matter affected 13 oscillations works perfectly well even in the limit Am 2 2 —> 0 , when solar effects become utterly 
unobservable {Toj. All procedures referred to above, however, rely fundamentally on a nonzero value of |C/ e 3 1 and 
neutrino oscillations in matter. In the next two sections, we argue that neither the former nor the latter is required 
in order to determine the mass hierarchy through neutrino oscillations — at least in principle. 


* The minimum value of | t/ e 3 1 2 required for unveiling the mass hierarchy via matter affected 13 oscillations depends on the specific 
properties of the neutrino source and the detector. For detailed discussions, see III 
' is the Vex. —► v /3 oscillation probability. In matter, antineutrino oscillation probabilities are given by the same expression as the 
neutrino ones, with A —> —A and U a i —*► I/*-. 

’ on 

+ This is required if we wish our two candidate hierarchies to fit all currently available atmospheric neutrino data. 
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To conclude this section, it is illustrative to mention that the “solar mass-hierarchy” was revealed by measuring how 
matter effects modify oscillations in the 12 sector. Since A m\ 2 is chosen positive-definite, the “solar mass-hierarchy” 
question is the following: is the v e flavor eigenstate “predominantly iq” or “predominantly 1 / 2 ,” be., is |[/ e i | 2 larger or 
smaller than | t/ e 2 1The electron-type neutrino survival probability in vacuum, assuming that Ato 2 3 23 effects are 
averaged out, is incapable of addressing this question. Under these circumstances, P ee is only a function of |f7 e i| 2 |Ue 2 | 2 
and hence invariant under |C/ e i | 2 «-> |t/e 2 | 2 . Matter effects break the degeneracy, such that, for 8 B solar neutrinos, 
P ee ~ | U e 2 \ 2 . Current data reveal that |C/ e 2 1 2 < 0.5 (and hence |f7 e i | 2 ~ 1 — |C/ e 2 1 2 > |C/ e 2 1 2 ), at more than the five 
sigma level. One reason we have been able to disentangle the solar mass hierarchy is that 9 12 is “different enough” 
from 0, 7 t/2 . Had that not been the case, we would have been unable to determine whether u e ~ v\ (9 12 = 0) or 
v e ~ V 2 (0i2 = 7t/2) — clearly a physical question! — via matter affected 12 oscillations. 

III. UNVEILING THE HIERARCHY THROUGH VACUUM OSCILLATIONS 

In the case of only two neutrino flavors, the neutrino mass hierarchy E3 is simply unobservable when it 
comes to studying vacuum flavor oscillations. This is easy to see. With only two flavors, P a p = 1 — P aa = 
sin 2 29 sin 2 {Am 2 L /4E), and there is simply no sensitivity to the sign of Am 2 or the sign of cos 29. This does 
not, of course, mean that the mass hierarchy is not physical — it is just not observable if one measures only vacuum 
oscillation probabilities. 

In the case of three or more flavors, this is no longer the case. Pure vacuum oscillations can, in general, determine 
whether the neutrino mass hierarchy, as defined in Sec. [fl is normal or inverted. In this section we explicitly present 
how this can be done in principle, and then review how it could be done in practice by comparing accelerator and 
reactor D e disappearance searches. 


A. How To Compare a Normal and an Inverted Mass-Hierarchy 


In vacuum, the survival probability P aa can be written as 

P aa = l-4|C/ al | 2 |C/ a 2 | 2 sin 2 -4|U al | 2 |t4 3 | 2 sin 2 - 4|I7 a 2 | 2 |f / a3 | 2 sin 2 , (H LI) 


where Ay = Ato 2 j -/ 2£’. Note that A 23 = A 13 — A 12 . Eq. lllll.ll) is valid, of course, independent of the mass hierarchy. 

We wish to probe whether normal and inverted mass hierarchies necessarily lead to different values of P aa . We 
will concentrate our efforts on comparing inverted and normal hierarchies which are capable of fitting the current 
neutrino data. More specifically, we will only consider normal and inverted mass hierarchies that share the same 
value of Am 2 2 and the same mixing matrix. These constraints clearly satisfy the requirement above and substantially 
simplify our presentation. Relaxing these constraints would render our analyses much more opaque, and would not 
modify qualitatively our results. 

Hence, we address the following question: for a given choice of Am 2 3 = Ato 2 ^ corresponding to a normal mass 
hierarchy (A rn 2 ^ > A m\ 2 > 0 ), is there an inverted mass spectrum characterized by the same mixing matrix and 
Am 2 2 but a different Am 2 3 = Am^ (—A m 2 ^ > Ato 2 2 > 0 ) that yields exactly the same survival probability P aa , 
for all values of L/El 

Define x by A m 2 ^ = — Ato 2 ^ +x. Note that x is restricted to be small enough so that A m 2 // is indeed a negative 
quantity.* In practice, x turns out to be of order Ato 2 2 , as we will see shortly. The difference between P aa (Am 2 g") 
and P aa (A T ri 2 3 ), P+ a - P~ a , is 


PL ~Paa = -4|U a3 | 2 \\U al 


\u a2 \ 2 


sm 


Ai 3 L 


— sin" 


. 2 ( (^13 — X)L 


2 ( (Ai 3 — A 12 )L^ . 2 ( (A 13 + A 12 — X)L^ 

sm ( - - - J — sm ( - - - j 


(HI-2) 


§ Since | C/ e 1 1 2 /1 1 2 = tan 2 6 12 , th is i s equivalent to inquiring about the sign of cos2$i2, or whether the solar angle is in the light or the 

dark-side of the parameter space M- 

* As a matter of fact, the magnitude of x should be small enough so that both Am^ and |Am^ 3 | agree with the currently allowed range 
of values for the “atmospheric” mass-squared difference. 
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where X = x/2E and, for simplicity, we have dropped the '+’ superscript in Amf^. There is, in general, no constant 
value of x that renders Eq. (inm exactly zero for all L.' The difference in the first square-bracket requires x = 0, 
while the difference in the second one calls for x = 2Am 2 2 . Curiously enough, there is a special solution in the case 
\U a i\ 2 = \U a2 \ 2 -. x = A mf 2 . 

Qualitatively, it is easy to understand what is happening. The distinction between a normal and an inverted mass 
hierarchy is, as far as vacuum oscillations are concerned, the following. For a normal mass hierarchy, Amf 3 > |Ato| 3 |. 
For an inverted one, the situation is reversed: |Aro 2 3 | < |Aro 2 3 |. One can choose, as is depicted in Fig.Q all normal 
and inverted mass hierarchy frequencies (|Ay|) to agree by equating |Ato 2 3 "| = lAm^l and |A? 7 i 2 3 | = |A? 7 i. 2 3 | 
equivalent to x = A m\ 2 - In this case, however, P+ a ^ P~ a , because the Ai 3 and A 2 3 amplitudes are, in general, 
different (except in the special case \U a i\ 2 = |f7 a2 | 2 )- 

We conclude that it should be possible to distinguish between the two hierarchies by studying neutrino oscillations 
in vacuum. The key requirement is the ability to tell A m 2 3 from A? 7 i | 3 or, in other words, simultaneous sensitivity 
to A m 2 2 and Am 2 3 effects. All experiments performed to date fail to have this simultaneous sensitivity for different 
reasons. “A m \ 2 experiments,” like KamLAND, lack sensitivity to Ato 2 3 effects, which “average out”.I “A m \ 3 
experiments,” like K2K, lack sensitivity to A m \ 2 effects, which are suppressed by a small A\ 2 L. 

From the discussion above, it seems that an experiment capable of measuring A m 2 3 with an uncertainty A(Am 2 3 ) < 
A m \ 2 should be able to tell a normal from an inverted mass hierarchy. However, it turns out that this is not 
automatically the case. 

Assume one is probing P aa with neutrino energies and baselines such that A 13 L ~ 1, as is the case for proposed 
next-generation accelerator and reactor neutrino experiments E3- Under these circumstances, Ai 2 P 1 and 

PL - P* a = —2L|/7 q3 | 2 sin (A 13 L) | [\U al \ 2 X + \U a2 \ 2 (-A 12 + (A - A 12 )] + O ((A 12 P) 2 , (XL) 2 ) . (IIL3) 


We are interested in values of x such that XL -C 1, and have made use of the fact that 

sin 2 (0 — e) = sin 2 9 — e sin 2(9 + 0(e 2 ). (III.4) 

One can easily find x such that, neglecting 0((Ai 2 L) 2 , (XL) 2 ) terms, P+ a — P~ a = 0: 


x = 


2|P q2 | 2 

|P a l| 2 + |Pa2| 2 


A?7l?2, 


(III.5) 


for |ZPq- 3| 2 ^ 0 and ignoring the possibility that sin(A 13 L) = 0. 

This means that if one probes, say, P w in vacuum at values of L and E such that A\ 2 L <C 1 and measures, 
assuming a normal mass-hierarchy, Ato 2 ^" = (2.13 ± 0.04) x 10 -3 eV 2 (a 2% measurement of the “atmospheric” 
mass-squared difference), there will be a statistically equivalent solution for an inverted mass hierarchy and Ara^ = 
— (2.02 ± 0.04) x 10 -3 eV 2 if a; MAI = 0.11 x 10 -3 eV 2 .§ Today, x^ is constrained to be 


Xfj./j, — 


WL 2 

1 - l^3l 2 


Am 2 2 


G 


[0.044,0.17] x 10“ 3 eV 2 , 


(III.6) 


taking three sigma upper and lower bounds for the oscillation parameters 

Appearance channels also contain a similar normal-inverted degeneracy. If one were to probe v a —> vp in a setup 
where A \ 2 L <C 1, a normal and inverted mass-hierarchy would yield the same oscillation probability P a p (up to 
0((Ai 2 i) 2 ) effects) if 


x a /3 = — 2cos£ 


U a 2 P/3 2 

U a 3U/33 


Ato? 2 , 


(HI-7) 


where ^ = arg[P a3 I/| 3 P* 2 P ; g 2 ]. 


t Keep in mind that x is restricted to be small, such that AmJj < 0. P + — P~ = 0 is guaranteed to be true for x = 2Am| 3 , independent 
of whether the oscillations occur in matter or in vacuum. This is trivial to understand, since x = 2Am^ 3 renders the two mass patterns 
in question identical —- such a large x value “inverts the inverted hierarchy.” 

* On top of that, \U e ^\ turns out to be too small. 

§ We define x aa , a = e,/q r as the value of x required to render the is a survival probability for an inverted and normal hierarchy 
indistinguishable in a range of values of L and E. Similarly, we define x a p as the equivalent quantity for Va. —> vp transitions 
(h,/ 3 = e,/x,r). 
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FIG. 2: Currently allowed values of x^/j, x x ee (top) and a:^ x x^t (bottom) (in units of A m\ 2 ) taking two (solid), three (dotted) 
and four (dashed) sigma allowed values of all oscillation parameters, according to iS). See text for details. The thin solid line 
corresponds to x e e,x^ T = 


B. Comparing P(u e —> v B ) to P(i/ M —> v^) 

One way to eliminate the degeneracy presented above is to compare two distinct oscillation channels, since the value 
of x that renders inverted and normal hierarchies indistinguishable depends on the parent and daughter neutrino 
flavors. For concreteness, we will discuss under what conditions one can determine the mass-hierarchy by studying P e 
and disappearance (at, say, a reactor neutrino experiment and a long-baseline conventional neutrino beam setup, 
respectively). One variant of this procedure has already been explored (using a different “language”) in the literature 
|13| and is, to the best of our knowledge, the only oscillation probe of the neutrino mass hierarchy identified so far 
that does not rely significantly on matter effects. In the case of a reactor neutrino experiment, the vacuum expression 
for the survival probability is a very good approximation. For accelerator-based experiments, this is only true for 
low-energy neutrinos (E < few hundred MeV) assuming that the neutrinos are propagating inside the Earth’s crust. 
Some consequences of nontrivial matter effects are discussed in Secs, [n] (for nonzero | C/ e 3 1 ) and II VI (for |C/ e3 | = 0). 

In the disappearance experiment, the two hierarchies can be made compatible by choosing x = x as defined 
in Eq. (IIII.6I) . On the other hand, for the reactor experiment, the two hierarchies will yield identical results if x = x ee , 
given by 


%ee 


2\U e2 \ 2 

1 - \U e3 \ 2 


Ami 2 


2Aro ? 2 sin 2 9 12 G [0.033,0.068] x 10 " 3 eV 2 , 


(III. 8 ) 


taking current three sigma upper and lower bounds for sin 2 $12 and A m\ 2 . As long as x w 7 ^ x ee , the “wrong hierarchy” 
measurements of A m 2 3 obtained via the electron-type antineutrino and muon-type neutrino disappearance analyses 
will disagree, thus breaking the degeneracy. Figureffltop) depicts the currently allowed values of and x ee , in units 
of Amf 2J if all oscillating parameters are allowed to vary within their two, three, and four sigma allowed ranges § 3 ]. 

In order to realistically determine the mass hierarchy, both the reactor and accelerator experiments are required 
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to measure the atmospheric mass-squared difference with great precision: A(Amf 3 ) < | x ^ — x ee \. Given what is 
currently known of the oscillation parameters, these could conspire in such a way that — x ee \ -C A 77 if 2 , as can 
be clearly seen in Fig. E^top). Improved knowledge of the oscillation parameters (especially 612 and $ 13 ) will help 
determine how challenging this proposition really is. 

While a 2% measurement of A rnf 3 in an accelerator u fJ disappearance experiment does not seem too far-fetched, 
as precise a measurement in a reactor experiment seems too optimistic. “Beta beams,” for example, may serve as a 
more advantageous alternative capable of precisely studying v e —> v e transitions in vacuum 0 . We are aware of an 
ongoing in-depth analysis of how to determine the neutrino mass hierarchy via a comparison of v tl —> v M and v e —> v e 
( v e —> v e ), to appear shortly El- 

Similar options which are potentially useful include comparing v fl —> u fl to ry, —> v T or —> v e . The currently 
allowed values of x fHl and x^ T are depicted in Fig. [SJbottom). Needless to say, precisely measuring A mf 3 in an 
electron or tau appearance channel seems far from realistic. Furthermore, as far as tau appearance experiments are 
concerned, neutrino energies are required to be significantly larger than 1 GeV, and the pure vacuum approximations 
we have been discussing so far simply do not apply in the case of an Earth-based accelerator experiment. We will 
briefly comment on the —> v T channel in the next section. 

Comparing the inverted-normal hierarchy degeneracy in —> and —> v e oscillations requires more care. The 
value of x^e obtained from Eq. (IIII. 71) diverges as \Uez \~ 1 as \U e z\ —► 0, and one may be tempted to conclude that 
x^e could end up much larger than rendering this comparison very sensitive to the mass hierarchy. This is not 
necessarily the correct conclusion for a few reasons. One is that, for small values of | C/ e 3 1, the A m\ 3 dependency of 
the ly, —> v e oscillation probability disappears, and at some point all sensitivity to “atmospheric” oscillations goes 
away — x^ e may be large, but not experimentally accessible! Another is that, in the limit x Amf 2 , Eq. Mil. 71 ) 
need not apply, since it was derived under the assumption that XL -C 1. 


IV. THE MASS HIERARCHY AND VANISHING \U e3 \ 

We emphasized in Sec. |H| that when |C/ e 3 1 ~ > 0 the traditional approach to determining the mass hierarchy through 
matter effects in the 13 sector fails, and new exploration channels are not only desirable but necessary. Furthermore, 
the approach studied in the previous section also fails in the limit |C/ e 3 1 ~» 0. One reason is simply that, in this case, 
oscillations involving electron neutrinos do not depend on Am 2 3 (including its sign).* Hence, in the limit |C/ e 3 1 —> 0, 
we are forced to explore the neutrino mass hierarchy through P MAl , P^ T and P TT . 

For vanishing $ 13 , neutrino oscillation experiments can still distinguish, in principle, an inverted mass hierarchy 
from a normal one. For example, in vacuum, P — P~^ (see Eq. cm) still does not vanish for all L/E for any choice 
of x if 0i3 = 0, unless $12 = 7 r/4 (maximal solar mixing, experimentally excluded). The same is true of P/j/ T — P// T - It 
is, however, curious that no new information is obtained if one compares the muon disappearance, tau disappearance, 
and tau appearance channels (for the same value of L and E). This can be seen in a variety of ways and boils down 
to the fact that, for U e 3 = 0, all three oscillation channels change “in the same way” when one compares a normal to 
an inverted mass hierarchy. For example, in the limit | t/ e 3 1 = 0 and A 12 A -C 1, 

Xfi[j, — Xt T Xf^T — 2A/71^2 COS 012, (-^^--0 

using Eqs. Mll.5Mll.7ll . The same behavior is observed in matter or outside the limit A 12 A <C 1. For details, see 
Appendix It suffices, therefore, to concentrate one’s efforts on the muon disappearance channel, which is what we 
do in this section. Throughout, we will assume |C/ e 3 1 = 0, unless otherwise noted. 


A. Comparing Values of P(;z M —> in Vacuum at Different L/E 

In Sec. I1I1I we argued that as long as experiments are performed under conditions such that A 12 L <C 1, for every 
normal hierarchy choice of mass-squared differences there is an inverted hierarchy one that yields the same oscillation 
probability (up to (A 12 A ) 2 corrections): Am^ = -Am^ + x , where x is given by either Eq. MII.5II or Eq. mo . 
In order to break this degeneracy, it is clear that one needs to explore setups where A 12 L ~ 1. 

Values of A 12 A ~ 1 imply very long baselines or very small neutrino energies. Given that we are considering 
experiments in which a muon is created, we need to steer clear from very low energies, and are forced to look at long 


One concrete manifestation of this can be found in Appendix^] 










FIG. 3: Values of that render P^{A\ 2 L/2) = P M(J (Ai2L/2) as a function of A 12 Z// 2 , for cos 2012 = 0.4 (sin 2 #12 = 

0.3). See text for details. The plot range is restricted to x ^ values such that |a; M(i | is of order A m\ 2 - 


baselines. Under these circumstances, it is quite likely that we will have to deal with relevant matter effects (short of 
considering experiments that shoot neutrino beams into space). Nonetheless, we first warm-up by discussing how one 
should go about disentangling the mass-hierarchy with muon disappearance in pure vacuum oscillations, and turn to 
matter effects in the next subsection. 

We start by making the following point: for a fixed L and E there are always values of x that renders P + — P~ = 0. 
Explicitly, P — P~^ = 0 translates (after a long series of trigonometric manipulations) into 


. / (2A 13 + X)L 

Sm l 2 


tan 


(A 12 — X)L\ 


) + (cos 2012 ) tan 


A12 L 


= 0 . 


(IV.2) 


The values of x that satisfy the trigonometric equation above are x 


— 2 Am 2 3 + 4wr ( 77 ), integer n, and 


x = A m\ 2 


4 E 

~T 


arctan ( cos 29\ 2 tan 



+ 4n7r 



(IV.3) 


for integer values of n and the appropriate definition of the inverse tangent function. The first set of solutions violates 
the restrictions we impose on the allowed values of x (as already pointed out, x = — 2Amf 3 transforms an inverted 
hierarchy into a normal one) and will be ignored. Values of x that satisfy Eq. itTOl) are depicted in Fig. [I] as a 
function of A\ 2 L/2, for cos26b 2 = 0.4. As expected (see Eq. HIV. Ill ), for Ai 2 L <C 1, x = 2cos 2 9i 2 Amf 2 for n = 0. 

Figure 01 reveals that, for a fixed L/E, there are values of x of order Aro 2 2 that render P + = P~ , even when 
A 12 L = 0(1). Indeed, for large enough values of A 12 T, there are several (discreet) values of x of order A m \ 2 that 
do the trick.t On the other hand, there is no constant value of x capable of solving P + = P~ for all L/E , as already 
extensively discussed. It is, therefore, clear that while a normal and an inverted mass hierarchy lead to distinct P M/1 , 
if —> v fl processes are studied around a fixed value of L/E , it is not possible to tell one mass hierarchy from the 
other. In order to break the degeneracy, more experimental information is required. 

For example, if one were to perform a very precise “short baseline” experiment characterized by small values of A 12 T, 
one would find two degenerate solutions for the “atmospheric” mass-squared difference: Am^ and —Am 2 3 " + a; s bi- 
On the other hand, a “long baseline” narrow-band-beam experiment with typical energy and baseline such that 
Ai 2 -b ^ 1 , would measure, assuming that the neutrino mass hierarchy is indeed normal, Arn 2 ^" and —A-t-aqbi- 
The mass hierarchy would be determined by combining both data sets as long as the precision A(Ato 2 3 ) with which 


t As A 12 L —> 0, all x values that obey Eq. Jl \ for n ^ 0 diverge, and only the n = 0 solution is acceptable. 
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both experiments measure Ato 2 3 is better than |x s bi — Xibi|- As an example, if A 12 L ~ 4 at the “long-baseline” site, 
l^sbi - Sibil ~ 0.8Amf 2 . 

Fig. ca can also be used to estimate how precise a broad-band-beam experiment in vacuum needs to be in order 
to establish the mass hierarchy. If one takes this experiment’s data and “bins” it in L/E , each L/E bin (or Ai 2 T 
bin) has to be capable of measuring A mf 3 with A(Amf 3 ) < Amf 2 - The same is true, of course, of experiments 
with variable baselines, most notably atmospheric neutrino experiments. We return to these issues and make a more 
realistic assessment of the situation in the next subsection. 


B. Comparing Values of P(—> v^) in Matter at Different L and E 


When \U e3 \ 
Numerically 


0, the importance of matter effects is estimated by comparing A 12 with the matter potential A. 


A 12 ( A m \ 2 \ / 100 MeV \ ( P \ 1 

A ' \ 8 x 10~ 5 eV 2 J \ E J \3g/cm 3 J 


(IV.4) 


We are interested in neutrino energies larger than 100 MeV, so that v^X —* yX' charged current processes are 
kinematically allowed, and typical Earth densities, which are around 2-5 g/cm 3 . 

As shown in Appendix^ P MM in matter can be written as P MM in vacuum with A^ and 0 12 replaced by effective 
quantities. Matter effects will modify neutrino oscillations (and our ability to determine the mass hierarchy) in two 
ways: 

• A 13 — A 23 = A (2 , so that the difference between the effective 13 and 23 oscillation frequencies will be different 
than Ai 2 . In the case of vacuum oscillations, this difference was recognized as the quantity that determines how 
precisely one needs to measure A 13 effects in order to establish the mass hierarchy; 


• sin20(0 = (Ai 2 /A^)sin26» 12 , so that the effective solar angle is different from its vacuum counterpart. In the 
case of vacuum oscillations, tan 2 612 measures the ratio of the amplitudes associated with 13 and 23 oscillations, 
and plays an important role in determining the value of x which renders P 3 ^ = P^ for a fixed L/E. 

Two limiting cases are particularly easy to describe. One is the limit of strong matter effects: |A| Ai 2 . Under 
these circumstances, A^ Ai 2 , and sin 2 2 9/% —► 0. Hence, 


pA/$> A12 
MM 


~ 1 — sin 2 2023 sin 2 


^ (A 13 — Ai 2 COS 2 012 )L ^ 


(IV.5) 


in which case x = 2 Amf 2 cos 2 0 i 2 renders P 3 ^ = P^. This can be understood in the following way. In the case of very 
strong matter effects and neutrino (antineutrino) oscillations, sin 2 9 ™ 2 —> 1 (cos 2 0 (^ —> 1 ), so that <-> ( v M <-> v^) 

oscillations are governed by only one oscillation frequency |l 6 l.§ Under theses circumstances, it is impossible to tell 
one mass hierarchy from the other. 

Another easy-to-describe limiting case is A” 2 L <C 1. Under these circumstances, we can use the results of Sec. EH 
to determine the value of x that renders P^ = without too much effort. We start by computing the auxiliary 
quantity y m , defined by A ^’" 1 = — A ^’ 7 " 1 + Y m ( Y m = y m /2E), which renders P MM (A^ 7 ™) = P M/i (A(" 3 ’ m ). The 


definition of y m is identical to the one for x (see Sec. IIIID after all parameters are replaced by effective ones. While 
y m has no simple physical interpretation, it can be related to x in a simple way and we can easily solve for it. 

Because the expression for P^ in matter is the same as the one in vacuum with the vacuum mixing parameters 
replaced by effective ones, we can simply read off the solution from Eq. 11111.511 : Y m = 2 cos 2 0 ^A 3 
A 12 cos 20 i 2 . We can relate y m to x: 


= A?*, -A- 


—A — A 12 + 2A 13 + A 


i 12 


A + A 12 - 2A+ - AT 


+ Y r 


A 13 = -A+a + A + A 12 - + Y r ‘ 

A 13 = A i3 + A', 


(IV. 6 ) 

(IV. 7) 
(IV. 8 ) 


^ A narrow A 12 L bin will contain several AmJ 3 “wavelengths,” which allows one to, in principle, perform a measurement of Am] 3 . We 
will return to this issue with more details in the next subsection. 

§ A13 — A12 cos 2 612 = A13 sin 2 0i2 + A 23 cos 2 0 i 2 - Hence, high energy (E > few GeV) v^ oscillations inside the Earth probe 

neither Am 2 3 nor Am 2 3 , but the solar-mixing weighted average of the two. This agrees with intuition obtained from the trivial limits 
012 —► 0 or 7 t/ 2, when oscillations can depend only on Am 2 3 or Am 2 3 , respectively. 
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FIG. 4: Values of x MM /Am 2 2 that render P^(L,E) = P^(L,E) in constant matter density p, as a function of L/E, for 
Am 2 2 = 7.9 x 10 -5 eV 2 , sin 2 #12 = 0.3, and p = 0 (vacuum oscillations, solid black line), L = 1000 km, p = 1.50g/cm 3 (solid 
red line), L = 3000 km, p = 3.14g/cm 3 (dashed pink line), and L = 6000 km, p = 3.73g/cm 3 (dotted blue line). See text for 
details. 


where in the last step we simply repeat the definition of x. Hence, 

x = ^ = Y m + A + A 12 -A? 2 , (IV-9) 

and, since Y m = — A + A 12 cos 29i 2 , x = 2A m\ 2 cos 2 6 * 12 . This result agrees with the vacuum value of x obtained 

in the case A 12 P <C 1, as expected. As is well known, for two flavor oscillations p™ tteT = P™ c in the limit A m L <C 1 
and it also turns out that here, in the limit A /\L <C 1, p™ tter = p™= for all a,(3 = e, p. r. 

The two limiting cases discussed above reveal that, in order to break the normal-inverted degeneracy in <-> u tJ 
oscillations, it is necessary to probe P MM for L and E values such that A\2 L > 1, independent of whether there are 
matter effects. Numerically, 


A12P 


/ Ato 2 2 \ /0.5 GeV\ / L \ 

\8 x 10~ 5 eV 2 J V E J V 3000 km / 


(IV.10) 


Furthermore, one needs to steer clear of very strong matter effects, which translates into E < few GeV for Earth- 
based sources and detectors and very long baselines. Similar neutrino energies and baselines have, of course, also been 
identified as the optimal conditions to study matter effects in the 12 sector m- 

We proceed to discuss how to determine the mass hierarchy following the strategy discussed in the previous subsec¬ 
tion. We compute the value of x that renders P+^ = P~^ at a fixed value of L/E. Following our previous discussion, 
it is straightforward to obtain 


AE 

x = 2 EA + AH-arctan 

12 L 


cos 2 055 tan 


A ? 2 L 


„ E . 

■^nn—, integer n, 
E 


(IV.ll) 


using Eqs. and lITWQll . 

Figure 0] depicts x values that render P^ = P^ in matter as a function of L/E for different baselines, assuming 
Antj 2 = 7.9 x 10 -5 eV 2 and sin 2 0\ 2 = 0.3. We choose matter densities equal to the average matter density value 
along the neutrino path, assuming that the neutrinos are produced and detected near the Earth’s surface, and using 
the preliminary Earth model [18j]. We extend all curves to very large values of L/E, even if these correspond to 
experimentally unexaminable values of the neutrino energy (E 100 MeV). In order to avoid an overly busy plot, we 
only display a few values of n (see Eq. (ITvTTli ) for each L and E. 

As advertised, significant effects are expected for large values of L (L > 1000 km) and small energy values (L/E > 
5000 km/GeV). In this region, not only is x ^ x(L/E —► 0), but also a; matter ^ x vac , as long as L/E is not too large. 
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Very large L/E values correspond to very small energy values and A 12 A, in which case matter effects become 
unobservable. This behavior can be seen in Fig.0] for L = 3000 km and L/E > several x 10 4 km/GeV. This regime 
cannot be observed in practice, as A 12 A implies E < 100 MeV for typical Earth matter densities. 

Fig.^Jreveals several distinct ways of determining the mass hierarchy. As briefly discussed in the previous subsection, 
a “broad-band” experiment with, say, a 6000 km baseline capable of precisely measuring the atmospheric mass- 
squared difference for energy values around, say, E = 1.2 GeV ( L/E = 0.5 x 10 4 km/GeV) and E = 300 MeV 
( L/E = 2 x 10 4 km/GeV) should be able to determine the mass hierarchy as long as A(Am 2 3 ) < 0 . 5 A? 7 i 2 2 . By 
“measuring the atmospheric mass-squared difference at some particular energy value,” we mean the following. Select 
a narrow window around a particular L/E value (say L/E±0.1L/E). Since we need values of L/E = 0( 10 4 km/GeV), 
we should find at least a few atmospheric oscillations inside our L/E window. Numerically, 


A i3^ = lf) 1 ( A?u 2 3 \ ( L ! E \ 

2 ^2.5 x IQ " 3 eV 2 ) \10 4 km/GeV ) ’ 


(IV.12) 


and it is clear that a very big challenge for “broad-band” neutrino experiments is to make sure that atmospheric 
effects do not “average out” — one needs to be able to “see” not the first and second oscillation maxima, but the 
tenth and the eleventh! This implies, among other things, excellent energy resolution (and angular resolution in the 
case of atmospheric neutrinos), and precise knowledge of the neutrino beam. A detailed quantitative study of the 
experimental requirements is clearly outside the scope and philosophy of this paper. 

One can also determine the mass hierarchy by performing two different experiments at the same value of L/E but 
different L. For example, a 3000 km baseline experiment performed at L/E around 2 x 10 4 would find x ~ 1.4Amf 2 
(among other values), while a different experiment performed around the same value of L/E but, say, L = 6000 km, 
would encounter x = 0.7A m\ 2 - If both experiments can determine the atmospheric mass-squared difference with an 
uncertainty smaller than 11.4 — 0.7|Am 2 2 , they should be able to unveil whether the neutrino mass hierarchy is normal 
or inverted. 


C. Comparing P(zq, —> u^) to P(z/ M —> v^) 

In the previous subsections, we described how different ty, disappearance searches can reveal the neutrino mass 
ordering, including the fact that neutrinos will often traverse a significant amount of matter on their way to the 
detector. When this is the case, it is well known that P/p will differ from P/p because matter effects violate CP- 
invariance^ (and CPT-invariance). It seems, therefore, plausible that one should be able to determine the mass 
hierarchy by comparing muon neutrino and muon antineutrino disappearance. 

We expect Ppp 7 ^ Ppp only for Ai 2 L > 1 and A ]2 /A ~ 1 given that (i) in the limit Ai 2 A or Ai 2 L <C 1, 
Ppp = P/u/ : • (ii) in the limit A Ai 2 , Ppp = Ppp, both given by Eq. (IIV.5I) . All of these features can be observed in 
Fig. GJtop), which depicts the values of x that render P 4 ^ = P and those that render P~ = Ppp as a function of 
L/E, for L = 2000 km, A m\ 2 = 7.9 x 10 -5 eV 2 , sin 2 9 12 = 0.3, and p = 2.83 g/cm 3 . 

Fig. Isjbottom) reveals that the comparison of results obtained with a muon neutrino and antineutrino beam 
performed at L = 2000 km and, say, E = 200 MeV (L/E = 1 x 10 4 km/GeV) will determine the neutrino mass- 
hierarchy as long as A(A m\ 3 ) < 0.5Ato 2 2 . 


V. SUMMARY AND CONCLUSIONS 

One of the challenges of next-generation neutrino experiments is to establish the neutrino mass hierarchy: is the 
relatively isolated z /3 state (much) heavier than zq and z/ 2 , or is it (much) lighter (see Fig. QJ? While we currently 
do not know the answer to this question, we do know that these two possibilities — a normal or an inverted mass 
hierarchy — are dramatically different.* In the case of a normal hierarchy, the neutrino masses can be ordered just 
like the quark and charged lepton masses — mi <C ?n 2 <C m3. In the case of an inverted mass hierarchy, the neutrino 
masses are ordered in a most unusual way — m3 <C mi , 2 and m 2 — m\ (m 2 + mi), i.e., two of the neutrino mass 


^ For [ 4,3 = 0, there are no intrinsic CP-invariance violating effects in neutrino<->neutrino oscillations. 

* We discount the extraordinary possibility that all three neutrino masses are quasi-degenerate. Whether the neutrino masses are all 
quasi-degenerate cannot be probed by oscillation experiments, and is already somewhat constrained by searches for neutrinoless double 
beta decay and indirect studies of the relic neutrino background. 
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L/E (km/GeV) x10 


FIG. 5: (top panel) Values of and x^/A m 2 2 that render, respectively, P/L(L,E) = P~^{L,E) (solid red line) and 

Pi(L, E) = Ppp(L, E) (dotted blue line), and (bottom panel) the minimum value of |a: MM — Xp.p\/A m 2 2 , as a function of L/E. 
L = 2000 km, p = 2.83g/cm 3 , A mf 2 = 7.9 x 1CF 5 eV 2 , and sin 2 #12 = 0.3. For comparison, the equivalent = Xpp, values in 
vacuum are also depicted in the top panel (solid black line). 


eigenstates are quasi-degenerate. We are yet to encounter two fundamental particles with quasi-degenerate masses. 
It is not surprising, therefore, that the determination of the neutrino mass hierarchy has been identified as one of the 
key discriminators of neutrino mass models ||]. 

In the literature, studies of matter-affected 13 oscillations have been singled out as the most realistic means of 
establishing the character of the neutrino mass hierarchy. It is interesting that these studies do not rely on properties 
unique to three (or more) neutrino flavor oscillations, and significantly mirror the (successful) path taken in order to 
determine the “solar” neutrino mass hierarchy. Such studies, briefly described in Sec. gh require \U e a\ to be large, 
and utterly fail if the limit | U e 3 1 —> 0 applies. 

In this paper, we have concentrated on probes of the neutrino mass hierarchy that do not rely on neutrino oscillations 
in matter or on \U e $\ being “large.” Instead, they rely heavily on the fact that, given that there are at least three 
neutrinos, P a p will depend on at least two different oscillation frequencies. 

First, we show that, in principle, vacuum oscillations are capable of distinguishing an inverted from a normal mass 
hierarchy — for a given normal mass hierarchy, there is no inverted mass hierarchy capable of producing the same 
P a p(L/E). We also argue that only neutrino oscillation experiments that are simultaneously sensitive to A mf 3 and 
Am'f 2 effects can “see” the mass hierarchy. Hence, we need experiments capable of measuring |A?n| 3 | with precision 
A(Ato 33 ) smaller than A m\ 2 - 

Second, we argue that experiments performed at L/E values such that A 12 L <C 1 suffer from a normal-inverted 
degeneracy, which can only be broken if more experimental information is available. The same is true of “narrow¬ 
band-beam” experiments which probe “around” a fixed value of L and E , regardless of the magnitude of A 12 T. 

If \U e 3 ,\ ^ 0, the neutrino mass hierarchy can be unveiled by combining searches for disappearance with searches 
for v e (or i/ e ) disappearance, as long as both experiments are capable of measuring Am'j 3 with precision A(Amf 3 ) 
smaller than A m\ 2 - It seems unlikely that this precise a measurement of A?nf 3 via v e disappearance can be performed 
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by reactor neutrino experiments. The KamLAND reactor experiment has measured Amf 2 with precision A(A?ti 2 2 ) — 
8%, which we take as an optimistic estimate of the precision with which a next generation reactor experiment can 
measure Amf 3 if |C/ e 3 1 2 > 10" 2 . 

In the limit \U e3 \ —> 0, any information regarding the neutrino mass hierarchy (if one sticks to neutrino oscillation 
probes) has to come from searches of or disappearance at distinct baselines and neutrino energies. If these 
experiments are to be performed on the Earth, at least one of them will require a very long baseline ( L/E > 
5000 km/GeV) and relatively small neutrino energies (E < 1 GeV) in order to guarantee that Ai 2 L > 1 and to 
avoid very strong matter effects ( \Z2GfN b A 12 ), respectively. Different potential scenarios include comparing 

disappearance searches at short (Ai 2 L -C 1) and long (Ai 2 L ~ 1) baselines, comparing searches for disappearance 
at different L but the same value of L/E , or comparing and disappearance in the presence of nontrivial matter 
effects (\Z2 GfN b ~ Ai 2 ). Note that “broad-band” accelerator experiments or studies of atmospheric neutrinos are 
potentially capable of performing all of these comparisons without the aid of other experimental setups — one can 
interpret a broad-band/atmospheric experiment as a collection of multiple experiments at different values of L/E. 

Instead of concentrating on the feasibility of any of the setups summarized above, we chose to identify the conditions 
that need to be met in order for one to successfully distinguish an inverted mass hierarchy from a normal one. We find 
that, irrespective of the scenario, we require multiple experiments capable of measuring |Ato 2 3 | at the few percent 
level (6.4% > Am 2 2 /|Am 2 3 | > 2.2%, taking three sigma bounds for mass-squared differences |(J). Whether this 
can be achieved in practice will be left for detailed experimental studies. We would, however, like to summarize 
the progress that can be expected in the foreseeable future.t Today, A(Am 2 3 ) ~ 25% ,6|. In the near future, the 
MINOS experiment expects to measure the atmospheric mass-squared difference with A(Am\ 3 ) ~ 10% [uj, while 
next-generation accelerator experiments aim at A(Am 2 3 ) ~ 1% |7]. In order to achieve such precise numbers, it is 
clear that improved knowledge of the neutrino beams and the neutrino-target cross sections are needed. 

We wish to conclude by reemphasizing that neutrino oscillation experiments can, in principle, address a most 
relevant and fundamental question in high energy physics — what is the character of the neutrino mass hierarchy? 
This capability does not depend on whether the oscillations occur in matter or on whether |C/ e 3 1 is nonzero. If, 
however, \U e3 \ is sufficiently small, our study here shows that unveiling the mass hierarchy through oscillations may 
prove very hard. Detailed studies are required in order to determine how hard. 
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APPENDIX A: ANALYTIC EXPRESSION FOR THE OSCILLATION PROBABILITIES IN CONSTANT 

MATTER DENSITY AND VANISHING \U e3 \ 


Simple analytic expressions for constant matter density neutrino oscillation probabilities can be obtained in the 
case \U e3 \ = 0. These are presented in this appendix. Other very useful approximate expressions, including 0{\U e3 \) 
corrections to the results derived below, can be found, for example, in 1201 . 

Neutrino oscillations in matter are described by the following set of differential equations: 

l ^Z = Ml 7 ’ ( A1 ) 

where 


v = 



(A.2) 


' Note that none of the proposed experiments currently under serious consideration meet the requirements for revealing the neutrino 
mass hierarchy spelled out above: either the baseline is too short, the energies too high, or the energy resolution insufficient to probe 
13 oscillations at low energies. 
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is a vector of neutrino flavor eigenstates and 

/ 0 0 0 \ / A 0 0\ 

M = U 0 Ai 2 0 U f + 0 0 0 . (A.3) 

V 0 0 a 13 j V 0 0 0 j 

Am 2 

The A ij = 2 ^ parameters contain all information regarding the neutrino mass-squared differences, and describe 

the energy dependence of the system. The forward scattering of electron neutrinos in matter is accounted for by 
A = V2G F N e , where G F is the Fermi coupling constant, and N e is the electron number density of the medium. As a 
function of the mass density p 1 A can be expressed numerically as A = 3.78 x 10~ 5 eV 2 g ^ m3 after we fix the electron 
to baryon fraction Y = 1/2. U is a unitary mixing matrix. In the standard parametrization [jj] with 0i3 set to zero, 
and ignoring potential Majorana phases, which are not relevant here, 

( cos 0 i 2 sin 0 i 2 0 \ 

- sin 0i2 cos 023 cos 012 cos 023 sin 023 = (ux u 2 U 3 ), (A.4) 

sin 0 i 2 sin 0 23 - cos0i 2 sin 0 2 3 cos 0 2 3 J 


where Ui are column vectors. Since U is unitary, u\uj = Sij, which means that Hi form an orthonormal basis for the 
space of three component vectors. 

The solutions to the three coupled differential equations in Eq. <ED can be expressed in terms of the eigenvalues 
A i and eigenvectors u) of M, i = 1,2, 3. Even though M is a three by three matrix, \ and vl may be found very 
easily. 

M, when U is given by Eq. iTPl) . may be expressed in the following convenient form: 

M = A 12 U 2 U 2 + A 13 M 3 U 3 + A(cos0i2Mi + sin0i2U2)(cos0i2Wi + sin0i2U2) t - (A.5) 


It is clear that one of the eigenvectors of M is M 3 : MM 3 = A 13 M 3 . Furthermore, since M is Hermitian, the other two 
eigenvectors are orthogonal to m 3 , which means they can be expressed as orthogonal linear combinations of m, and m 2 . 
The problem of finding the three eigenvalues and eigenvectors of a three by three matrix is now reduced to finding 
two eigenvectors in a two by two effective matrix, which is very easily solved. Indeed, the problem reduces to the well 
known case of two-flavor neutrino oscillations in constant matter, and its solution is very well known 

Hence, the resulting oscillation probabilities are conveniently expressed in terms of effective mixing angles and 
frequencies, which are simply related to their vacuum counterparts. 


m 

12 — 

\JA? — 2AAi2 cos 20±2 H~ ^12 ? 

(A. 6 ) 

m 

—A — A 12 + 2 A 13 + A ^2 

(A.7) 

13 

2 

m 

- A - A 12 + 2 A 13 - A% 

(A. 8 ) 

‘23 

2 

I m 

12 ~ 

Ai 2 . „„ 

7vT sm20 i2> 

^12 

(A.9) 

im 

A 12 cos 2012 — A 

(A.10) 

12 

A- 

S3 

II 

023, 

(A.ll) 


A^ — A^g = A^, and A”^ reduces to A 12 in the vacuum limit A 12 A. 

The oscillation probabilities coincide exactly with the corresponding vacuum expressions once the vacuum oscillation 
parameters are replaced by the corresponding effective ones. We explicitly list all oscillation probabilities in constant 
matter below, keeping in mind that, due to T-invariance,* P a p = Pp a , a, (3 = The oscillation probabilities of 

antineutrinos are given by the same expressions, with A —A. 


P p „ = 1 — sim 2671 sin 2 


LAV 


(A.12) 


When | C/ e 3 1 = 0, there are no T-invariance violating effects in neutrino oscillations in vacuum. Furthermore, a constant matter potential 
does not induce T-invariance violating effects in neutrino oscillations in matter. Matter-induced CP-invariance violation is, of course, 
present. 
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(A.13) 


(A. 14) 


As expected, P ea (a = e, fi, r) does not depend on Amf 3 , and hence does not care whether the neutrino mass hierarchy 
is normal or inverted. This is not the case of which are given by 



(A.15) 


(A.16) 


(A.17) 


The terms of P Wi(1TiTT that depend on Amf 3 are identical, up to an overall sign. Hence, |P + — P~\ is the same for 
all non-electron oscillation probabilities. Any normal-inverted degeneracy observed in a study of —> v tJ will be 
present in exactly the same form in a similar study of v tJ —> v T or v T —> v T . 
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